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We study a spin system with both two- and four-spin 
exchange interactions on the triangular lattice as a possible 
model for the nuclear magnetism of solid ^He layers adsorbed 
on grafoil. The ground state is analyzed by the use of the 
mean-field approximation. It is shown that the four-sublattice 
state is favored by introduction of the four-spin exchange in- 
teraction. A possible phase transition at a finite temperature 
into a phase with the scalar chirality is predicted. Applica- 
tion of a magnetic field is shown to cause a variety of phase 
transitions. 
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I. INTRODUCTION 

Nuclear magnetism of the two-dimensional ^He layers 
adsorbed on a grafoil has been a puzzle since Greywall 
and Busch reported the. specific heat data at the cov- 
erage p = 0.178 A~^.El'El At this coverage the adsorbed 
■^He atoms form two layers and each layer solidifies into 
a triangular lattice. The first layer is a high-density solid 
and is considered not to contribute to the magnetism 
in the mK temperature region. The specific heat data 
down to 2 mK showed a prominent peak at 2.5 mK as is 
expected from the conventional Heisenberg antiferromag- 
net (HAP) on the triangular lattice. But integration of 
their data extrapolated linearly to lower temperatures in- 
dicated that a half of the total magnetic entropy will not 
be released down to the absolute zero temperature. Elser 
proposed that the second layer forms a \/7 x sjtructure 
in registry with the first layer triangular lattice.El Then it 
decomposes into inequivalent A and B sites which may 
be magnetically decoupled. Since A sites form a kagome 
lattice, he proposed that the missing entropy may be ex- 
plained by the HAF on the kagome lattice rather than 
that on the triangular lattice. Inspired by this proposal 
many investigations were accomplished on the HAF on 
the kagome lattice in recent years.u According to numer- 
ical studies of finite clusters the kagome HAP shows a 



second peak in the specific heat at a low temperature, 
which might compensate the missing entropy.Lli3 On the 
other hand, Roger showed for a cluster with 16 spins 
that a second peak can be also produced by effects of the 
multi-spin interactions which-, should be of considerable 
amount in solid "^He systems J3 

Recent specific heat data down to O.lmK at the cov- 
erage p — 0.180 A~^, however, do not exhibit any sharp 
second peak iruthe temperature region where the kagome 
HAF predicts. Q Furthermore the estimate of the entropy 
released at temperatures above O.lmK nearly exhausts 
the total magnetic entropy of the second layer. The peak 
at T ~ 2.5 mK found by Greywall and Busch turned 
out to have a much broader tail in the low temperature 
region than was first thought. It is clear that these fea- 
tures of recently obtained specific heat data cannot be 
explained in terms of either the HAF on the triangular 
lattice with nearest neighbor exchange couplings or that 
on the kagome lattice. So the behavior of the specific heat 
remains as a puzzle though there is not missing entropy 
any longer. 

It is well known that magnetic interactions due to ring 
exchanges involving more than two spins are important 
in solid '^He. The famous uudd state in bcc solid ^Jie is 
incorporated by the four-spin exchange interaction.El The 
exchange frequencies in bcc solid '^He computed by using 
the path-integral Monte Carlo me^od agree fairly well 
with estimates from experiments. □ For a monolayer of 
■^He at the coverage p = 0.0785 A~^ on grafoil the cyclic- 
exchange frequencies were computed by Bernu et al. in 
a similar way.E3 Their result shows that the three-spin 
exchange term predominates the others at this coverage 
and the high-temperature specific heat is governed by the 
two-, three- and four-spin exchange interactions while the 
six-spin exchange term as well is necessary to be taken 
into account to estimate the Weiss temperature properly. 

The recent experimental results indicate that we have 
to take account of all spins in the second layer in order 
to understand its magnetism. We assume that the dif- 
ference between A and B sites causes negligible effects 
on the translational symmetry of the magnetic interac- 
tions and that the magnetism of the second layer may 
be described by an = 1/2 spin system on the regu- 
lar triangular lattice. We employ the Hamiltonian with 
the two-, three- and four-spin exchange interactions as 
an effective Hamiltonian for the second layer of adsorbed 
■^He, though no estimates of exchange frequencies were 
reported for the second layer so far. This model may de- 
scribe the '^He monolayer adsorbed on plated graphite as 
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well.tllll3 Since three-spin exchange terms map into usual 
two-spin exchange coupUngs, the model is expressed by 
the two-parameter Hamiltonian 

n.n. 

H^jY, T,-crj+KY,hp, (1) 

<i,j> p 

where cri/2 is the nuclear spin of the '^He atom at the site 
i and denotes the sum over all possible diamonds 
composed of two nearest neighbor unit triangles. The 
four-spin exchange term of the diamond p depicted in 
Fig. 1(a) reads 

hp— ^ (Ta ■ (Tp + {(Ta ■ (Tb){(Tc ■ (Td) 

<Q,/3> 

+ {(^a ■ (^d){o-b ■ O-c) - {a-a ■ (Tc){(Tb ' CTrf), (2) 

where X]<a /3> denotes the sum over all possible pairs of 
four spins on the diamond. The parameter J and K are 
related to the exchange frequency J„ of the n-body cycle 
of the nearest neighbors as 

J - - J2/2 + J3 (3) 

and 

K = - J4/4. (4) 

It should be noted that tiie exchange frequency J„ is 
always negative for any The magnitude of the two- 
body exchange frequency J2 is expected to be large at a 
low enough atomic density of the Iffffr but it decreases 
with density more rapidly than Ja.Q'tl Thus the param- 
eter J is expected to be positive at a low density while 
it should be negative in the high density region. On the 
other hand K is always positive. We expect that a wide 
parameter region of the ratio of J to if may be realized 
by varying the coverage of ^He layers in a realistic sys- 
tem. It is therefore useful to investigate what kind of 
the ground state occurs in the system described by the 
Hamiltonian (|l|). Calculations of the specific heat and 
the susceptibility of a Ifi-spin cluster were reported by 
Bernu et al. previously J13 

II. THE GROUND STATE IN THE CLASSICAL 
LIMIT 

As a first step to study the spin system described by 
the Hamiltonian we examine the ground state in the 
classical limit. Since we neglect quantum fluctuations, we 
treat cTj as a unit classical vector in the following. When 
K = Q the ground state of the system is well known, 
that is, the ferromagnetic state for J < and the state 
with the so-called 120° structure for J > 0. On the other 
hand we cannot guess the ground state of the system 
intuitively when J = 0, since the Hamiltonian is governed 
by the unusual four-spin term. So we use a variational 



argument. We divide the total Hamiltonian into a sum 
of those of interpenetrating hexagons with 7 spins. That 
is 

h 

where Hh is the Hamiltonian of the hexagon depicted 
in Fig. 1(b) and is expressed as a sum of Hamilto- 
nians (H) for six diamonds contained in the hexagon. 
Then we searched numerically for the spin configuration 
which minimizes Hh. It has turned out that the lowest- 
energy spin configuration is a four-sublattice spin struc- 
ture where any two spins on different sublattices make 
the relative angle a where cos a = —1/3 (see Fig. 2(b)). 
We can divide the whole lattice into four sublattices such 
that four vertices of every diamond belong to different 
sublattices. The spin configuration which minimizes a 
hexagon, therefore, can be extended to the whole lattice 
in such a way that the state has the lowest energy for 
all hexagons in eq. (||) consistently and so realizes the 
ground state of the whole lattice. The spin configuration 
has the four-sublattice structure, where the spin vectors 
point the vertices of a tetrahedron if we put their bot- 
toms at its center, as shown in Fig. 2(b). Hence we call 
this spin structure as "tetrahedral" structure. In fact we 
can prove in the same way that the tetrahedral struc- 
ture is the ground state at least in the parameter region 
-K/2 < J <2K. 

For general values of J and K we determine the ground 
state in the mean-field approximation considering both 
three- and four-sublattice spin structures. In a three- 
sublattice structure with the sublattices a, b, and c, the 
energy is expressed as 

E/N = ( J + 4:K){(Ta -(Tb + fTb-fTc + a-c- (To) 

+ 2K{{(Ta ■ (Tb){(Tc ■ (Ta) + {(Tb ■ Crc){(Ta ■ (Tb) 

+ {(T,-(Ta){(Tb-(Tc)} + iK. (6) 

In a four-sublattice structure with the sublattices a, &, c 
and d, we have 

E/N ='^{J + (iK)Y,f^o.-'yii 

a<l3 

+ K{{cra ■ (Tb){(Tc ■ (Td) + {cTa ' Crc)((Tb • CT d) 

+ (cTa ■ (Td){crb ■ CTc)}. (7) 

In above equations ctq denotes the spin vectors on the 
sublattice a etc. We determined the ground state spin 
configuration by comparing the energies of the two con- 
figurations which minimize eqs. (||) and (|^), respectively. 
Minimization of eqs. (^ and (^ were done numerically. 
As a result we obtained four ground state phases with 
following spin structures: 

1. the 120° structure with E/N = -3(J + K)/2 for 
Q<K < 3J/25, 
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2. the tetrahedral structure with E/N = -J-lTK/i 
for < 3J/25 < X and < -3J/8 < K, 

3. the uuud structure with E/N = —3K for < 
-J/8 <K < -3K/8, 

4. the perfect ferromagnetism with E/N — 3J + 21K 
for < X < - J/8. 

The uuud structure has four sublattices and spins on the 
three sublattices are ahgned parallel while spins on the 
other one is antiparallel to them. The ground-state spin 
structures are depicted in Fig. 2. 

The most interesting and novel one among the ob- 
tained ground-state spin structures is the tetrahedral 
structure. It is remarkable that this structure has a scalar 
chiral long-range order. The scalar chirality in the unit 
cell may be defined as 

XS ^O-a- {a-b X (Tc) + (Tb ■ {(Td X (T c) 

+ (Tc ■ {o-d X (Ta) + (Td ' (cTf, X Cr„), (8) 

which is an order parameter of Ising-type with a discrete 
Z2 symmetry. Since the tetrahedral ground state has a 
finite value of xsi 16V3/9, the system may undergo a 
phase transition of Ising universality at a finite tempera- 
ture in spite of the isotropic rotational symmetry inherent 
in the Hamiltonian. This contrasts with the 120° struc- 
ture which has a vector chiral order and may cause piily 
topological phase transitions at finite temperatures. t£l 

III. PHASE DIAGRAM IN THE MAGNETIC 
FIELD 

Application of a magnetic field to frustrated spin sys- 
tems is known to cause interesting magnetization pro- 
cesses and/or various phase transitions. First we con- 
sider the case with K — and J > 0, i.e. the HAF on 
the triangular lattice with nearest neighbor coupling J. 
It is well known that two phase transitions occur under 
the magnetic field at T=O.EJ When the magnitude of the 
field H is less than 3 J (the magnetic moment of a spin 
is assumed as unity), The ground state is the so-called 
umbrella state with a coplanar spin configuration where 
the spins on a sublattice orient downwards (the magnetic 
field is assumed to be applied upwards) and those on the 
other two sublattices orient upwards but canted to the 
field by an angle (see Fig. 3(a)). The angle decreases 
with the magnetic field and vanishes at iJ = 3J. For 
3J < H < 9J, spins on all sublattices are canted to 
the field but two of them are parallel to each other (see 
Fig. 3(c)). We call the phase with this spin structure as 
C3 phase. Spins are perfectly polarized for H > 9J. At 
H = 2>J the spins are collinear, i.e. those on two sub- 
lattices are oriented upwards while those on the other 
sublattice point downwards (see Fig. 3(b)). We call this 
spin structure as uud structure. It was shown that the 
phase with the uud structure expands and is realized in 



a finite range of the magnetic field aided by the thermal 
or the quantum mechanical fluctuationsO'c3 

We searched for the ground state for general values of 
J and K numerically minimizing eqs. (|6[) and (^ added 
by -[al + al + al)H/i and -(a- +af +af + a^)i//4, re- 
spcctively. We found seven ground-state phases in total. 
They are the following: 

1. the umbrella phase, which appears for a weak mag- 
netic field when J > and K is small. 

2. the uud phase which extends around H ~ 3 J(> 0) 
and for small K. 

3. the C3 phase, which appears for a wide range of the 
parameters both for positive and negative J near 
the boundary of the fully polarized phase. 

4. the weak field four-sublattice phase, the spin struc- 
ture is an analogue of the umbrella structure but 
is not coplanar. The spins on a sublattice point 
downwards and those on the other three sublattices 
orient upwards but canted to the field by an an- 
gle (see Fig. 3(d)). There is a three- fold rotational 
symmetry about the magnetic field axis. We call 
this structure "mushroom" structure. This phase 
appears for large K and small H. 

5. the uuud phase which appears for small | J| and 
rather weak field. 

6. the high field four-sublattice phase with the struc- 
ture analogous to the C3 phase. Spins on three sub- 
lattices are parallel to each other and all spins are 
canted to the field (see Fig. 3(e)). We call this 
phase as C4 phase. 

7. the ferromagnetic phase where the spins are per- 
fectly polarized. This phase appears for H > 
9{J + 8K). 

Above spin structures are shown in Fig. 3. The phase 
diagram is depicted in Fig. 4 for J > and in Fig. 5 for 
J < 0. We adopted a scale as J -I- 8K = 1 for J > and 
— 2 J + 8K = 1 for J < for convenience. The full lines 
denote the first order transitions while dotted lines the 
second order ones. The transitions between the three- 
and the four-sublattice structures are necessarily of the 
first order. The transitions between the structures with 
the same sublattice structures are of the second order in 
most cases. An exception is the boundary between the 
mushroom and the uuud phases for J < 0, where the 
phase change occurs with a magnetization jump. The 
critical field He at the second order transition points are 
determined analytically and tabulated in table 1. 

There is a tetracritical point at He ^ QJ ^ 24K when 
J = AK. Magnetization process may show peculiar be- 
haviors around this point. As an example magnetization 
per spin m is shown in Fig. 6 as a function of the applied 
field for J = 0.4 and K = 0.75. In this case the sys- 
tem undergoes six phase transitions with variation of H. 
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There appear two plateaus in the magnetization curve at 
m = 1/3 and m = 1/2. The 1/3- and 1/2-plateaus cor- 
respond to the uud and uuud phases, respectively. Fur- 
thermore the curve shows magnetization jumps at the 
first-order phase-transition points. Another magnetiza- 
tion curve is shown in Fig. 7 for the parameter J — —0.1 
and K = 0.1. 



IV. SUMMARY AND DISCUSSION 

We have shown that various ground-state phases oc- 
cur on the triangular lattice by coexistence of the two- 
and the four-spin exchange interactions. We have found 
new phases with four sublattices, i.e. the tetrahedral 
structure and the uuud structure. Since we have as- 
sumed that the ground states have either the three- or 
the four-sublattice structure, there remains a possibility 
that a state with a longer period may be more stable 
than those obtained above. In some cases, however, we 
are certain that the present results are rigorous in the 
classical limit. That is, the case when K = (the pure 
Heisenberg model), and the case when J = and H = 0. 
So it is quite reasonable that the obtained phases are sta- 
ble in the regions close to these limiting cases. In fact we 
could prove that the tetrahedral phase is stable for small 

It is remarkable that the phase with scalar chirality ap- 
pears when the four-spin exchange dominates. This sug- 
gests the existence of a phase transition of the Ising uni- 
versality at a finite temperature. Our preliminary Monte 
Carlo simulations in the classical limit clearly exhibit di- 
vergence of the specific heat at T ~ X for J = 0. The 
present model seems to be the first example of a real- 
istic spin system in two dimensions with full rotational 
symmetry which undergoes a chiral phase transition. De- 
tailed analysis of th&iinite-temperature properties will be 
reported elsewhere.ll3 

In the region where the competition between the two- 
and the four-spin interactions are serious, more study is 
necessary to be conclusive. In fact we found in Monte 
Carlo simulations that various states with long periods 
and zero magnetization have quite close energy to that of 
the uuud structure. Furthermore, it can be shown that 
a new 120° structure state with nine sublattices has the 
same energy as the uuud structure. This indicates that 
the frustration of the system is strong in the phase with 
uuud structure due to the competition between the two- 
and the four-spin interactions. The application of the 
magnetic field would select the uuud structure from the 
degenerate ground states, since it has the largest magne- 
tization among them. We hence believe that the ground- 
state phase diagram under the magnetic field (Fig. 5) 
will not be changed seriously by the appearance of long- 
period states. 

Effects of quantum fluctuations, which we have com- 
pletely neglected in the present study, are surely very 



important since the magnitude of the nuclear spin is 1 /2 
for ^He. It is still not quite established that even the 
well-known 120° structure is stable against quantum fluc- 
tuations when K — and J > (HAF), though recent 
numerical studies seem to support for the stabilityJl3 It 
is possible that the coexistence of the two- and the four- 
spin interactions destabilizes the ordered ground state 
and helps the quantum fluctuations to destroy it. The 
stability of the ground states against the zero-point spin 
wave fluctuations are preseiitly under investigation and 
will be reported elsewhereJlj 

Finally we compare our results with the existing ex- 
perimental ones, though it is still not quite clear whether 
the model is relevant to describe ^He layers. The specific 
heat data at low temperatures exist only at two differ- 
ent ■^Qilujes of the coverage, i.e. p = 0.180 and 0.228 
A~^.ElE3 No sharp peak was observed at both coverages. 
The system with p = 0.228 clearly corresponds to a 
negative J since it shows about 60 percepi of full polar- 
ization under a very small magnetic field.tll On the other 
hand the parameter J of the system with /9 = 0.180A~^ 
is still unclear. At this coverage the system has been 
thought to be antiferromagnetic since it has a negative 
Weiss temperature 0. However a negative 6 does not 
imply a positive J since 6 = — 6( J + 6K) holds in the 
present model. If the system corresponds to a positive 
J, the ground state at this coverage should be the pos- 
itive large-J phase or a new quantum phase caused by 
competitions between the two- and the four-spin intearc- 
tions and quantum fluctuations, since the system does 
not show a finite-temperature phase transition. If this is 
the case, there may occur the phase with a chiral long- 
range order between two coverages p = 0.180 A~^ and 
0.228 A~^ (if the density of the second layer varies con- 
tinuously), and we may observe a sharp peak in the spe- 
cific heat by varying the coverage of the adsorbed layers. 
Another possibility is that the system with p — 0.180 
A~^ has a negative J whose value corresponds to the 
uuud ground state in the mean field theory. In this pa- 
rameter region we have seen that the strong frustration 
leads to the existence of various states with long peri- 
ods whose energies are quite close to that of the ground 
state. Then the effects of quantum fluctuations may re- 
alize a novel ground state without magnetic long-rage 
order. Observations under a magnetic field will be use- 
ful to determine the phase that the system belongs to. 
Above discussion is based on the assumption that more- 
than-four-spin exchange interactions are not important. 
If this is not the case it is necessary to investigate a spin 
model with more-tlian-four-spin exchange interactions in 
search for possible ground states. 
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Phasel 



Phase2 



He 



umbrella 
uud 

ferromagnetic 

mushroom 

uuud 



uud 

C3 
C3 

uuud 

C4 



3(J - AK) 
3( J + 4i^) 
9( J + 8K) 

4:Ij + 2K) 

4:{J + 8K) 



TABLE L He at the second order transition point. 



Figure Captions 

FIG. 1. (a) A unit diamond described by the four-spin ex- 
change Hamiltonian (2). (b) A unit hexagon described by H^ 
in eq. (5). 



FIG. 2. The ground state spin structures: (a) the 120° , (b) 
the tetrahedral and (c) the uuud structure. The ferromagnetic 
spin structure is an obvious one. 



FIG. 3. The ground state spin structures in the magnetic 

field: (a) the umbrella, (b) the uud, (c) the C3, (d) the mush- 
room and (e) the C4 structure. The uuud structure is shown in 
Fig. 2. The magnetic field is assumed to be applied upwards. 



FIG. 4. The ground state phase diagram for J > 0. The 
parameters J and K are scaled as J -I- 8 A' = 1. The phases 
are labelled by u(umbrella), m(mushroom) and f (perfect fer- 
romagnetism) . Other labels are self-evident. The full and the 
dotted lines indicate the phase boundaries of the first and the 
second orders, respectively. 



FIG. 5. The ground state phase diagram for J < 0. The 

parameters J and K are scaled as —2 J + 8K = 1. The labels 
and the lines are the same as in Fig. 4. 



FIG. 6. The magnetization curve for J = 0.4 and K = 0.75. 

The system undergoes three first-order transitions and three 
second-order ones. The C3 phase appears twice. 



FIG. 7. The magnetization curve for J = —0.1 and 
K = 0.1. The system undergoes two first-order transitions 
and two second-order ones. 
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